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Abstract 

In this paper,we present a rigorous demonstration and discussion of the quantum adiabatic theorem for 
systems having a non degenerate continuous spectrum. A new strategy is initiated by defining a kind of 
gap, "a virtual gap", for the continuous spectrum through the notion of eigendifferential (Weyl's packet) 
and using the differential projector operator. Finally we obtain the validity condition of the adiabatic ap- 
proximation. 

PACS: 03.65.Ca, 03.65.Ta 

The adiabatic theorem is one of the basic results in quantum theory [TJ [2] . It is concerned with quantum 
systems described by an explicitly, but slowly, time-dependent Hamiltonian. There has been a sudden regain 
of interest in the adiabatic theorem for itself among physicists when in 1984 M . V. Berry [3] pointed out that 
if it was applied to Hamiltonians satisfying H (ti) = H (£2), it could generate a phase factor having non trivial 
geometrical meaning. And more recently, the adiabatic theorem has renewed its importance in the context of 
quantum control [1], for example, concerning adiabatic passage between atomic energy levels, as well as for 
adiabatic quantum computation [5]. 

There are several points of view for a discussion of the quantum adiabatic theorem; each one offers interesting 
insight. As T. Kato [2] has pointed out, the contents of the adiabatic theorem embody two parts: first, the 
existence of a virtual change of the system which may be called an adiabatic transformation, and, second, the 
dynamical transformation of the system goes over to an adiabatic transformation in the limit when the change of 
the Hamiltonian is infinitely slow. The adiabatic theorem proof given by M. Born and V. Fock [lj, although very 
general, is still restricted by the assumption of considering the purely discrete and non-degenerate Hamiltonian's 
spectrum, except for accidental degeneracy caused by crossing. These limitations are rather artificial from the 
physical point of view and should be removed from Kato's derivation of the adiabatic theorem. Several authors 
[5] had formally extended the Kato's results on the approximate validity of the adiabatic theorem when the time 
T, during which the approximation takes place, is large but finite. G. Nenciu [7J demonstrated the adiabatic 
theorem for bounded Hamiltonians. Later, J. E. Avron and A. Elgart [8] showed that the adiabatic theorem 
holds for unbounded Hamiltonians as well and applied it to deal with the quantum Hall effect. 

Let us simply recall here that the works following that of Born and Fock [T] by Kato [2] , Garrido [5] , Nenciu 
[7J and J. E. Avron et al. [7J have led to a formulation of the adiabatic theorem under the usual gap assumption 
9nm (t) — E n (t) — E rn (t), between level n and m. One may then state that a general validity condition for 
adiabatic behavior is well controlled as follows: the larger is the quantity < n jfi 1 [ffnm (t)] the smaller will be 

the transition probability. 

Despite the existence of extensive literature on rigorous proofs of estimates needed to justify the adiabatic 
approximation [2 [7J O E] , doubts have been raised about its validity [TU] leading to confusion about the precise 
condition needed to use it [11]. In part, this is because some papers emphasize different aspects, such as 
the asymptotic expansion, the replacement of the requirement of non-degenerate ground state by a spectral 
projection separated from the rest of the spectrum, dependence of first order estimates on the spectral gap, and 
even extensions to systems without a gap. Adiabatic theorem without gap conditions is know to be true [5], 
however, in general, no estimates on the error terms are available. J. E. Avron and A. Elgart have shown in 
ref. [8] that the adiabatic theorem holds provided the spectral projection is of finite rank independently of any 
spectral considerations. A similar result was proven by F. Bornemann [12J for discrete Hamiltonian when the 
set of eigenvalues crossings is of measure zero in time. The limitation of these approaches is that, in general, no 
estimate can be made on the rate at which the adiabatic regime is attained [5] . The gap condition is generally 
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associated to spectral stabilities. Consequently, the situation where the gap does not exist will led to spectral 
instabilities. Thus it is difficult to establish smooth spectral projections which is a necessary condition for 
the validity of the adiabatic theorem in the practical applications. In fact, the generalized adiabatic theorem, 
according to J. E. Avron and A. Elgart's approach [8J, is much more appropriate for the systems without a gap 
condition and which have a discrete origin. 

In this letter, we present a straightforward, yet rigorous, proof of the adiabatic theorem and adiabatic 
approximation for systems whose Hamiltonian has a completely continuous spectrum supposed non-degenerated 
for reasons of simplicity and which checks a certain number of conditions which will be given later on. 

In the case of continuous spectrum we cannot numerate eigenvalues and eigenfunctions, they are characterised 
by the value of the physical quantity in the corresponding state. Althoug the eigenfunctions ip (k; t) of the 
operators with continuous spectra cannot be normalised in the usual manner as is done for the functions of 
discrete spectra, one can construct with the ip (k; t) new quantities - theWeyl's eigendifferentials (wave packets)- 
[T5] which possess the properties of the eigenfunction of discrete spectrum. The eigendifferentials are defined 
by the equation 

k+Sk 

IMM)> = J \<p(k';t))dk'. (1) 

k 

They divide up the continuous spectrum of the eigenvalues into finite but sufficiently small discrete regions of 
size 8k (see FigfTJ). 



Figure 1: Decomposition of the continuous spectrum. 

The eigendiffercntial |T]) is a special wave packet which has only a finite extension in space; hence, it vanishes 
at infinity and therefore can be seen in analogy to bound states. Furthermore, because the Sip have finite spatial 
extension, they can be normalized. Then in the limit 5k — > 0, a meaningful normalization of the function ip 
themselves follows: the normalization on S functions. 

For 5k, a small connected range of value of the parameter k (this corresponds to a group of "neighboring" 
states, see FigfT]), the operator 

k+8k 

5P{k;t)= J \ip(k';t))(ip{k';t)\dk' (2) 

k 

represents the projector (the differential projection operator [14] ) onto those states contained in the interval 
and characterized by the values of the parameter k within the range of values 5k. The action of 5P (k; t) on a 
wavefunction \ip (t)} is defined by 



k+Sk 

5P{k;t)\Tp(t)) = J C{k'-t)\ip{k'-t))dk'. (3) 

k 

The application of the differential projection operator 5P (fc; t) causes thus the projection of the wavefunction 
onto the domain of states ip (fc; t) which is characterized by k values within the 5k interval .Before proceeding 
further, we give the statement of the adiabatic theorem. 
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Figure 2: Evolution of a range of energy of width 8k as a function of time. 



Let us call Ut (s) the evolution operator where s is the fictitious time and T is the time interval during 
which the evolution of the system takes place 

ih^-U T {s)=TH{s)U T {s), (4) 
us 

and the slowly time-dependent Hamiltonian H (s) = J E(k,s) \<p(k,s)) {(p(k,s)\dk , < s < 1, has a purely 
continuous spectrum E (k,s). 

If the following conditions are fulfilled 

(i) As it is mentioned earlier (see FigJT]) the continuous spectrum is divided into discrete regions of size 8k, we 

must define or create a gap of energy for the continuous spectrum, in other words, the size 8k is chosen 
so that 

E (k; s) - E (k'\ s) » i , Vfc' g [k, k + Sk] . (5) 

(ii) We assume that the eigenvalues are piecewise differentiable in the parameter s, and there is no level crossing 

throughout the transition (see Figj2j), in other words: 

E(k';s) ^ E(k";s) /a € [0, 1] , k' 6 [k, k + 8k] , k" £ [k, k + 8k] . (6) 

(iii) The derivatives -§^8P (k; s) and -§^zSP (k; s) are well defined and continuous in the interval < s < 1. 

Under these conditions it is possible to prove the adiabatic theorem: 

Theorem 1 If the quantum system with time- dependent Hamiltonian having a non degenerate continuous spec- 
trum is initially in an eigenstate \ip (k, 0)} of H (0) and if H (s) evolves slowly enough then the state of the system 
at any time s will remain in the interval [k, k + 8k]. 

The adiabatic theorem can be formally written, at the first order, in terms of the evolution operator as 

Vfc: lim U (s) 8P (k; 0) = 8P (k; s) lim U (s) + O ( )- J . (7) 

Notice that if, initially, the system is in the state \tp(k,0)) so that H (0) \(p(k, 0)) = E{k, 0) \<p(k, 0)) and 
expanding an arbitrary state vector on the basis of the instantaneous quasi-eigenfunction, then (J7|) implies 

lim U (s) \tp (fc, 0)} = 8P (fc; s) lim U (s) \ip(k; 0)) (8) 

T — >oo T — >oc 

k+Sk 

and in the limit T — > oo the state U (s) \<p (fc; 0)} = J C (fc'; s) \ip (fc'; s)) dk 1 , belongs to the subspace generated 

k 

by the states \tp (fc; s)) pertaining to the interval [fc, fc + 8k] . 
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Proof. The demonstration that we present follows the same approach developed in ref. [14] for the discrete 
case. To this effect, we introduce a unitary operator A (s) having the property 



SP (k, s) = A (a) SP (k, 0) A+ (s) Vfc G 5ft. 
It is completely defined by the initial condition A (0) = I and the differential equation 

d 



(9) 



(10) 

(11) 
(12) 

(13) 

The unitary transformation A + (s), applied to the operators and the vectors of the Schrodinger's picture, 
produces a new picture: the rotating axis picture: 

(14) 
(15) 



ih—A(s) =K(s)A(s). 
os 

The operator K (s) obeys the following commutation relation: 

ih^- s 5P{k,s) = [K(a),SP(k,a)], 
and is determined without ambiguity if we add the following supplementary condition: 

(ip(k;t)\K(t)\(p(k';t)) = 0, Vfc' G [k,k + Sk], 
equation that yields the following expression 

K (t) =ih [ [1-5P (fc; t)] \<p (fc; t)) (tp (fc; t)\ dk. 



(s) = A+ (s) H (s) A{s) = J E (fc, s) \<p (fc, 0)) (tp (fc, 0)\dk, 

similarly K ^ (s) becomes 

K {A) (a) = A+ (s) K (s) A (s) . 
The evolution operator in this new "representation" is U^ A ' (s) = A (s) Ut (s) -It is defined by 



ihlLuW (s) = \THW (s) - (a)] (a) , (0) = /. 

os L J 



(16) 



Since (s) and (s) are T-independent, it is to be expected that in the T — * oo limit the first term of 
the right hand side in (|16p dominates. We can approximately solve such an equation in which we go over to a 
time-dependent reference frame following the axis which diagonalize (s). We define $t (s) via 



ih^T (s) = THW (s) $ t (s) 

OS 



in which the solution may be written, with the initial condition $^ (0) = /, as 

iTa (fc, s) 



$ T (s) 



exp 



\ip(k,0))(ip(k,0)\dk, 



(17) 



(18) 



where a (fc, s) = J Q S E (fc, a') ds' . If, as we will see immediately, (s) tends toward $t (s) for large T, we will 
have approximately 



U T (s) 



A(s) $ T (s). 



(19) 



We go over to a second picture and we show that the remaining evolution operator differs from the identity by 
terms O (i) . Thus, we change to a last picture with operator!^ (s) = ^ (s) A + (s) Ut (s), and (— ) generator 
K (s) = $+ (s) A+ (a) K (a) A (s) $ T (s): 



d 

ih—W (a) = K (a)W (a) , W(0)=I, 
os 



equivalent to the integral equation 



W{s) 



ft Jo 



K{s')W (s')ds'. 



(20) 



(21) 
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Now, we prove that in the limit T^oo,W(s) = I + (4) .We begin by considering the operator F (s) 
J*q K (s 1 ) ds' . Any operator (and in particular F (s)) admits the following decomposition 



F(s) 

Using (fl8|) we obtain: 

F(k,k',s) 



exp 



F (k, k! , s) dkdk' 

(</?(fc,0) \K(s')\cp(k',0)) \<p(k,0)) ((p(k',0)\dkdk'ds'. 



K {A) {k,k's')ds' k'<£[k,k + 5k] 



iT (a(k,s')-a(k',s')) 



h 



(22) 



(23) 



an expression in which we have introduced the condition k' ^ [k,k + Sk] because, from ([12"]) . we deduce 
F (k, k', s) = for k' e[k,k + 8k] . 

Let k! (fi [k, k + Sk], since .K"^ (fc, fe's') is a continuous function of s, our assumption implies that a (k, s') — 
a (k' , s') is a continuous nonvanishing monotonic function of s; after integrating (123j) by parts we obtain 



F(k,k',s) 



h 



exp 



exp 



iT (a (fc, s') - a (k 1 , s')) 



(k,k's') 



iT (a (k,s')-a(k',s')) 



E(k,s')-E(k',s') 
E(k,s')-E(k',s') 



ds' 



(24) 



hence, according to the condition ([5]), F (k,k' , s), fc' ^ [fc, A; + <SA;] , asymptotically converges toward as y. 
Summarizing, as T — > oo we have: 



(25) 



Using (|2U|) . integration by parts turns (f2"Tj) into: 

W (s) = / + l -F (s) W (s) + \ 



F(s')K(s)W(s')ds', 



(26) 



since the last two terms in this equation contain the factor F(s), then for T — > oo and from Ut (s 
A (s) 4> T (s) W (s) we obtain 

U T (s) =s A (s) $ t ( s ) L " /J 



T 



(27) 



Finally (HU) implies <J> T (s) <5P (fc, 0) = SP (k, 0) $ T (s) and hence A (s) $ T (s) <5P (k, 0) = A (s) SP (k, 0) $ T (s) = 
SP (k, s) A (s) <3>t (s). This concludes the proof of the adiabatic theorem ([7]). ■ 

If T is sufficiently large, we can, in first approximation, replace U (t±,to) by its asymptotic form: 

U (h,t ) = U T (I) ~ A (l)<f> T (l). (28) 

This is called the adiabatic approximation. If the initial normalized state is \ip (ko, 0)), under this approximation 
U(ti,to) \ip(ko,0)) ~ A(1)$t (1) |y(^0)0))- To determine the validity of the adiabatic approximation for a 
given process, we can estimate the error by computing the probability rj of finding the system at time t\ in a 
state different from A (1) $ T (1) \<p (k , 0)}: 



rj = (feo, 0) |f/ + (ti, * ) 4 (1) $ T (1) QoA + (1) $+ (1) U (t ls io) | p (fe , 0)) 
where Qo = I — SP (ko, 0). This quantity may be rewritten as 

= (ko, 0) |TU+ (1) Q 1U (1)| <p (k , 0)) . 
Solving (|21[) iteratively and keeping only the first order term, we find 

V « ^2 <^(fco,0)|F+(l)Q o P (l)|^(fc ,0)) 



(29) 
(30) 



fc0[fca,fco+5feo] 



|(^(fc ,0)|P(l)|^(fc,0))| 2 &. 



(31) 
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Now, let us define a normalized time through the variable transformation t = to + sT (0 < s < 1), and the 
initial normalized state \ip (ko, to)) of H (to) with the eigenvalue E (ko, to)- Then, using (fT3"|) , (j2"3"|) and performing 
the change s — > t in eq. (j3"Tj) yields 



(^(*b,*o)|F(ti)| ¥ j(fc,t )> = ifi 



ex Pi - 



[£(fco,i')-£(M')]^' 



in 



[< V Oo,t) |£ (fc, t)) - (*o, t) |<SP (*d,t)| (p (k, t))] dt. 
Since (ip(ko,t) \SP(ko,t) \ (p(k,t)) = for k [k Q ,k + Sk ] equation (PTj) may be recast as 



(32) 



1 



ITi 



exp 



in 



[E (k ,f) - E (k, t')\ dt' \ (<p (k , t) \<p (k, t)) dt 



to 



' kf[ko,ko+Sko] 

and the adiabatic approximation for \ip (kg, to)) holds only if r\ <C 1 which requires 



dk, 



(33) 



k.L) 



iTi 



In 



ex Pi - 



[E (k , t') - E (k, t')] dt' \ {ip (k 0l t) \ip (k, t)) dt 



<1, Vfc ^ [k ,k + Sk ] ■ 



(34) 

the integral of eq. (|34|) will be sufficiently small, if the phase of the integrated function vibrates fast enough 
and the amplitude of the integrated function is small enough, thus ^p(fc ^fc ! t) is at the maximum 



Sp( ko ^k,t) » max 
iei*0i*ij 



h(ip(ko,t)\(p(k,t)) 



E(ko,t') -E(k,t') 



, Vfc ^ [ko, k + Sk ] 



(35) 



The condition 77 <g; 1 is therefore, in most cases, certainly satisfied if max ^p(fe ->fe t) *C 1 or equiva- 

k^[ko,ko+Sk ] 

lently, 

max \((p(ko,t)\<p(k,t))\<. min \E (k ,f) - E (k,t')\ , Vte[to,h], (36) 

k£[ka,k +Sk ] k^[k a ,k +Sk ] 

with max and min taken over all k [ko, ko + Ska]. Condition (|36[) may be taken as a criterion for the validity of 
the adiabatic approximation in the case of a continuous spectrum. This estimate of the adiabatic approximation 
could not be made in the Avron-Elgart's approach [S] as mentioned earlier. 
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